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The nonequilibrium dynamics of carriers and phonons in graphene is investigated by solving 
the microscopic kinetic equations with the carrier-phonon and carrier-carrier Coulomb scatterings 
explicitly included. The Fermi distribution of hot carriers are found to be established within 100 fs 
and the temperatures of electrons in the conduction and valence bands are very close to each other, 
even when the excitation density and the equilibrium density are comparable, thanks to the strong 
inter-band Coulomb scattering. Moreover, the temporal evolutions of the differential transmission 
obtained from our calculations agree with the experiments by Wang et al. [Appl. Phys. Lett. 96, 
081917 (2010)] and Hale et al. [Phys. Rev. B 83, 121404 (2011)] very well, with two distinct 
differential transmission relaxations presented. We show that the fast relaxation is due to the rapid 
carrier-phonon thermalization and the slow one is mainly because of the slow decay of hot phonons. 
In addition, it is found that the temperatures of the hot phonons in different branches are different 
and the temperature of hot carriers can be even lower than that of the hottest phonons. Finally, we 
show that the slow relaxation rate exhibits a mild valley in the excitation density dependence and 
is linearly dependent on the probe-photon energy. 



PACS numbers: 78.67.Wj, 63.20.kd, 71.10.-w, 78.47.J- 

I. INTRODUCTION 

The unique electrical and optical properties of 
graphene, e.g., large mobility, long coherent length and 
exceptionally low electrical noise, make it a promising 
material for the development of nanoscale devices The 
performance of many such devices depends critically on 
the dynamic properties of carriers and phonons. There- 
fore, a thorough understanding of these properties is es- 
sential. 

The time-resolved optical pump-probe measurement 
is a powerful tool widely used to probe the ultrafast 
dynamics of photoexcited carriers and has been ap- 
plied extensively to graphene lately^— In these works, 
a fast differential transmission (DT) relaxation of sev- 
eral hundred femtoseconds, followed by a slower picosec- 
ond relaxation were observed**"— Dawlaty et al& sug- 
gested that the fast one is due to equilibrating of carriers 
through the carrier-carrier scattering and the slow one 
is related to the cooling of the hot-carrier distribution 
through the carrier-phonon scattering. However, both 
experimenta l 10 ! 11 and theoretical worksi^ show that the 
photoexcited carriers lose most energy to optical phonons 
within the time scale of the fast relaxation (about 500 fs) 
and the slow relaxation rate is in the same order of the 
hot-phonon decay rate obtained from the time-resolved 
Raman spectroscopy! 13 ! 14 Therefore, the fast relaxation is 
supposed to be associated with the rapid carrier-phonon 
thermalization and the slow one to the relaxation of hot 
phonons through the phonon-phonon scattering. 7-9 How- 
ever, the theoretical investigations on this problem in 
the literature are performed by using the coupled rate 
equations which calculate the energy transferred among 
the carriers, phonons and environment 2& This method 



is based on the ansatz that the carrier-carrier Coulomb 
scattering is very strong so that the carrier Fermi dis- 
tribution can be established very rapidly. This should 
be verified especially at high photo-excited carrier den- 
sity in which the screening is strong. Moreover, in their 
model, the temperatures of electrons in the conduction 
and valence bands are assumed to be always identical 
after buildup of the Fermi distribution. This is reason- 
able when the excitation density is much larger than the 
equilibrium carrier density, since the distributions of elec- 
tron and hole are almost identical in this case. Neverthe- 
less, it should be examined when these two densities are 
comparable. In addition, in the previous investigations, 
the contribution of the remote-interfacial (RI) phonons 
is neglected and the carrier-phonon scattering matrices 
for the longitudinal and transverse optical phonons near 
the r point are set to be identical*^— The influence of 
the missing physics, especially the different temperatures 
in various phonon branches from different carrier-phonon 
scattering strengths, is still unclear. All these questions 
suggest that a detailed theoretical investigation from a 
microscopic approach is essential. 

In this paper, we investigate the nonequilibrium dy- 
namics of carriers and phonons in graphene via the mi- 
croscopic kinetic equation approach with the carrier- 
phonon and carrier-carrier Coulomb scatterings explicitly 
included. The temporal evolutions of the carrier distri- 
bution and the phonon number as well as the DT are 
obtained numerically. We find that the hot-carrier Fermi 
distribution is established within less than 100 fs. Fur- 
thermore, due to the strong inter-band Coulomb scatter- 
ing, the temperatures of electrons in conduction and va- 
lence bands are shown to be very close to each other even 
when the excitation and equilibrium densities are com- 
parable. It is also shown that the calculated DTs have 
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good agreement with the experimental data by Hale et 
alZ and Wang et al£- for different graphene layer num- 
bers and excitation densities. Moreover, our calculations 
provide strong evidence to the claim in the previous ex- 
perimental works 7- — that the fast relaxation of the DT 
is due to the carrier-phonon thermalization and the slow 
one mainly comes from the hot-phonon decay. Further- 
more, it is shown that due to the different carrier-phonon 
scattering strengths, the temperatures of hot phonons in 
different branches are different. Moreover, the tempera- 
ture of carriers can be even lower than that of the hottest 
phonon. Finally, it is discovered that the slow relaxation 
rate exhibits a mild valley in the excitation-density de- 
pendence and depends linearly on the probe-photon en- 
ergy. 

This paper is organized as follows. In Sec. II, we set up 
the model and lay out the kinetic equations. In Sec. Ill 
the results obtained numerically from the kinetic equa- 
tions are presented. We summarize in Sec. IV. 



II. MODEL AND FORMALISM 

We start our investigation from graphene on Si02 or 
SiC substrates. Exploiting the nonequilibrium Green's 
function approach, the kinetic equations of the carriers 
can be constructed as 18-21 



dtffj,ki> — dtffiku\ cc + dtf^kv\ cp + <9t/ Mk „| 



(1) 



Here /i = +(— ) represents the K(K') valley, k stand for 
the wave vectors relative to the K or K' points and f^v 
represent the electron distribution functions in the con- 
duction {y — +) or valence {y = — ) bands. 9t/ M ki/| ee de- 
scribe the carrier-carrier Coulomb scattering terms and 
<9t/ A tki/| ep give the carrier-phonon scattering terms in- 
cluding the scatterings between carriers and acoustic, op- 
tical as well as RI phonons. In this investigation, we as- 
sume that the initial carrier distribution after the pump- 
ing is isotropic. Therefore, the carrier-impurity scatter- 
ing terms dtf^k^l ■ are always zero. The carrier-carrier 
scattering terms can be written as (ft = 1 throughout 
this paper) 

9tfiJ,Uu\ cc = — 47T ^ ^ -fky.k'y'-fki+k-k'i^.kii^ 
kV/V kii^J^ 

x l^kX' | 2< 5( £ k'i/' — £ki/ + £kj+k-kv; — £ ki^) 

X ( f > f < f > f < — f < 

/i'ki+k— V!v' x J /i'kii/ 2 lite- v ■> /iki; J /Vki+k— kVJ 
X //^ki^/^k'y'/^ki/)' (2) 

in which J k ^ kv = \[1 + vv' cos(6» k - 9 k ,)}, f< ku = f^ v , 
fuku = 1 — /pki/ and £ki/ = wpk with v-p being the Fermi 
velocity. V^ k , denotes the screened Coulomb potential 
under the random phase approximation 2 ^— 

K+^m = 0[i - ek +q- - £k -')]= (3) 



in which V® — 2-KVpr s jq is the two-dimensional bare 
Coulomb potential with r s being the dimensionless 
Wigner-Seitz radius 2 ^— and n(q, w) is given b y 24 ' 28 " — 

n(q, w )= 2W +q *' / T" / ^-n+ - (4) 

\iw ' k 

The carrier-phonon scattering terms are given by 



X (//VkV'//jk!v n k-k'A //^'k'i/'/^ki/ n k-k'A) • 



(5) 



Here A is the phonon branch index and w q ,\ is the cor- 
responding phonon energy; ri^r x = rt q A + \ ± § with n q A 
representing the phonon number. For acoustic phonons, 
w qAC = ^phQ with v p h being the acoustic phonon velocity 
and the scattering matrices are 
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-^kiy.k+qi/'<5, 



(6) 



in which D is the deformation potential and p m denotes 
the graphene mass densit y 31 ' 32 For the RI phonons, 



|M RI « 1 ' 1 2 



-2qd 



-^klAk+qi/' &fj 



aq 



|1- V q °n(q,£ kly -£ k4 



(7) 

where g represents the dimensionless coupling parameter 
depending on the material of the substrate ; 27 ' 33 a is the 
C-C bond distance, d stands for the effective distance 
of the substrate to the graphene sheet . 23 ' 25 " — For the 
optical phonons, 
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In this investigation, we include the transverse optical 
phonons (Kto) near the K(K') point and the longitudinal 
(Flo) as well as transverse optical (Fto) phonons near 
the r point. The corresponding parameters read 



^ vlkk , = (D 2 K )[l-^'co S (e k - 



ek'-2^_ k )]V.(9) 
M]<^,-/*' (10) 



with k = 1(-1) for r LO (r TO ) phononsi 5 -^ 

In this paper, the dynamics of the RI and optical 
phonons are studied, while the acoustic phonons are 
always set to be at the environment temperature To- 
We further adopt the assumption following the previous 
worki 7 -^ the phonons in the same branch can equilibrate 
themselves very quickly. Similar to those of the carriers, 
one can obtain the kinetic equations of the hot phonons 



d t nqx = <9 t 7i qA | op + d t riqx 



ipp 



(11) 



Here, dtriq\\ cp and <9t?i q A| pp come from the carrier- 
phonon scattering and the anharmonic decay of hot 
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phonons, respectively. For RI and optical phonons, they 
are given by 

ph i / / i 
^ kq' vv fifj, 

x l^+ q v'l 2 [/^/^k +q ^("q'A + 1) 

— //^ki>//I'k+qV n q'A] i (12) 



TABLE I: Parameters used in the computation. 



<9 t n qA | pp = 



J qA - %A 



(13) 



In above equations, t pp is the phenomenological relax- 
ation time from the phonon-phonon scattering; n qA is the 
number of the A branch phonons at environment temper- 
ature To; iVph = m(E max /vF) 2 /4:n is the number of the 
phonon modes participating in the carrier-phonon scat- 
tering, where E m&K represents the upper energy of the hot 
carriers which are able to emit phonons£&21 and m = 1 
for the RI, T T o and T L o phonons and m = 2 for the two 
degenerate Kto phonons at the K and K' point. 

By numerically solving the kinetic equations [Eqs. ([1]) 
and (|TT1) 1 with the same numerical scheme laid out 
in Ref. l35l the temporal evolutions of the carrier dis- 
tribution and the phonon number can be obtained. 
Then, the evolution of the optical transmission at 
the probe-photon energy uj pl can be calculated from 

T pi (u pT ) = |1 + A r iayC r (^pr)v / Mo/eo/(l + ^rcf)!" 2 where 
n rc { is the refractive index of the substrate and iVi ay 
is the number of graphene layers. The optical conduc- 
tivity is given by cr(w pr ) = -e 2 (f fJ ,u. ul + - / Mklo _)/4 with 
|ko,| = Upr/2v F MSrM The DT is then calculated from 
AT pr /T p ° r = (T pr - T p ° r )/T p ° r , with T p ° r representing the 
transmission before pumping. It is noted that the epi- 
taxial multilayer graphene can also be described by our 
model, since they have been demonstrated to have simi- 
lar phononic and electronic properties to those of single- 
layer graphene ^ 40 ' 41 Nevertheless, the differences be- 
tween these two systems, i.e., the carrier- RI phonon scat- 
tering becomes negligible when the number of layers is 
large, should be taken into account. The material pa- 
rameters used in our calculations are listed in Table Q] 



III. RESULTS 

In this section we first study the buildup of the hot- 
carrier Fermi distribution in Sec. IIIA. We show that 
the Fermi distribution with identical temperature in the 
conduction and valence bands can be established within 
100 fs. Then, in Sec. IIIB we simply use the hot-carrier 
Fermi distribution as the initial carrier distribution and 
compare the calculated DT with the experimental data. 
The evolutions of carrier and phonon temperatures are 
also investigated here. Finally, we study the excitation- 
density and probe-photon-energy dependences of the 
slow DT relaxation rates in Sec. IIIC. 
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FIG. 1: (Color online) (a) vlog(f~ ku — 1) as function of \eia,\ 
for conduction-band electrons [y = +) at t = 3, 30 and 100 fs 
and valence band electrons (y = — ) at t = 100 fs. The results 
at t — and 3 fs for conduction-band electrons are plotted in 
the inset, (b) Temporal evolution of temperatures of electrons 
in conduction (blue dots) and valence bands (red squares) 
from the calculations with (solid curves) and without (dashed 
curves) the inter-band Coulomb scattering filter- 



A. Buildup of hot-carrier Fermi distribution 



We set the initial carrier distribution to be 



(14) 
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in which F{e ku ) = 1/{1 + exp[(e k , - fj, ) / (k B T )]} 
is the carrier Fermi distribution before pumping with 
fx denoting the initial chemical potential; G(e-k v ) = 
Aexp [— (|eio/| — w pu /2) 2 /(2S 2 )] is the photo-generated 
carrier distribution with A and S representing the ampli- 
tude and standard deviation, respectively and uj pu stand- 
ing for the pump-photon energy. In our computation 
here, the substrate is chosen to be Si02 and the equilib- 
rium carrier density Nq is taken to be 6 x 10 11 cm~ 2 . We 
also set A = 0.15, S = 19.3 meV and w pu = 1.5 eV, cor- 
responding to the excitation density N ex = 8 x 10 11 cm -2 
and the absorbed intensity I a = 1.9 mJ/m 2 . Then the 
conduction-band-electron density is about 1.75 times of 
the hole density (hole distribution is defined as /^kh = 
1 — f/ik-)- The other parameters are taken to be To — 
300 K, £ max = 0.9 eV and Tpp — 3.8 ps. We first focus 
on the distribution of the conduction-band electrons and 
plot the evolution of log(f~£ + — 1) in Fig. [TJa) . Note that 
if the distribution / M ki/ is the Fermi distribution, one has 

log(/^ - !) = - ^)/(k B T v ), (15) 

with T„ and \i v representing the temperature and chemi- 
cal potential in the corresponding band. Therefore, if the 
curves in the figure become linear with \e\av\-, the buildup 
of the Fermi distribution is identified. In the inset of 
Fig.QJa), one finds a valley located at \s\J\ = 0.8 meV in 
the initial distribution, coming from the photo-generated 
carriers. With the evolution of time, the valley is rapidly 
smeared out by the carrier-carrier Coulomb scattering as 
shown by the curves with t = 3 and 30 fs in Fig. HJa). 
Moreover, the curve with t = 100 fs becomes almost lin- 
ear with |£ki/|, indicating the buildup of the Fermi dis- 
tribution. It is noted that this time scale is in the same 
order as those in the experiments in graphite. 1 • We stress 
that the hot carrier Fermi distribution presented here is 
obtained directly from the microscopic kinetic equations 
and no ansatz is needed as all the scatterings are in- 
cluded. 

Then we turn to the valence-band electrons. Our calcu- 
lations show that the Fermi distribution of the valence- 
band electrons is established in the same time scale as 
that of the conduction-band ones, and we only plot 
-log(/^_ - 1) at t = 100 fs in Fig. QJa). More im- 
portantly, one finds that the slopes for conduction (black 
dotted curve) and valence (yellow dashed curve) bands 
are very close to each other at t = 100 fs. This indicates 
that their corresponding temperatures T + and T_ are 
very close to each other [see Eq. (|15p], even though the 
conduction-band-electron density is about 1.75 times of 
the hole density. To make this more pronounced, we plot 
the evolution of the hot-carrier temperatures T u fitted 
from Eq. (fT5|) in Fig. [TJb). The difference between T + 
and T_ is shown to be less than 10%. This phenomenon is 
due to the strong inter-band Coulomb scattering, which 
can be seen by comparing the temperatures from the cal- 
culations with (solid curves) and without (dashed curves) 
the inter-band Coulomb scattering Hf* tei . 



B. Temporal evolutions of DT and temperatures of 
carriers and phonons 

As shown in the previous subsection, the hot-carrier 
Fermi distribution is established very rapidly and the 
temperatures of electrons in the conduction and valence 
bands are almost identical. Therefore, in the following 
calculations, the initial carrier distribution is set to be 

Wi = 0) = 1/{1 + exp[( £k , - £)/{k B T° e )]}, (16) 

where T® denotes the hot-carrier temperature; [P v repre- 
sent the chemical potentials in conduction [y = +) and 
valence [y = — ) bands. Tj? and can be determined by 
the equilibrium carrier density and the excitation density 
N ex as well as the absorbed intensity I a . For simplicity, 
we set I a = N ex uj pu with u; pu denoting the pump-photon 
energy. With this initial carrier distribution, the evolu- 
tion of the DT can be obtained numerically. 

We first compare the DT from our calculations with 
the experimental results by Hale el al£ in single-layer 
graphene on SiC>2 substrates [Fig. [U(a)].— Here T Q = 
300 K, w pr = 1.1 eV and u; pu = 1.5 eV as indicated 
in the experiment. In this case, the excitation density 
is much larger than the equilibrium carrier density, thus 
the chemical potential before the pumping has little in- 
fluence on the evolution of the DT and is set to be at the 
Dirac point for simplicity. Then Tj} = 4163 K and f/\_ = 
-fj,°_ = -478 meV correspond to N ex = 4.6 x 10 12 cm~ 2 
and I a = 11 mJ/m 2 which are close to the estimated 
values given in Ref. |7|. The fitting parameters here are 
£ mM = 0.9 eV and t pp — 3.8 ps. Our results agree very 
well with the experimental data and show a fast relax- 
ation with the characteristic time about 0.28 ps, followed 
by a slow one with 1.33 ps. 

To reveal the underlying physics of these two relax- 
ations, we plot the evolution of carrier and phonon tem- 
peratures in Fig. [2Kb). The carrier temperature T e can 
be obtained by fitting log(/^ CJ/ — 1). The temperature 
of hot phonons in A branch can be obtained from T\ = 
wa/[&b hi(l + l/riqx)]. From Fig. [D^b), it is seen that 
the temperatures of the phonons first increase rapidly 
and then decrease slowly, with the peaks very close to 
the crossover point between the fast and slow DT re- 
laxations [see Fig. Ufa)]. Since the fast increase of the 
phonon temperatures is due to the rapid equilibration of 
the carrier-phonon system (less than 500 fs)^2r— and the 
decrease comes from the slow hot-phonon decay (about 
several picosecond) ) 13 ' 14 the fast and slow relaxations of 
the DT can also be attributed to these two processes, 
respectively^ This result supports the conjectures in the 
previous experimental works In addition, by com- 
paring the slow relaxation of the DT with the exponen- 
tial fitting curve [black dotted curve in Fig. Ufa)], one 
finds that the relaxation rate slightly increases with the 
temporal evolution when t > 4 ps. This can be under- 
stood via the approximate formula Eq. (IA.4[) , which can 
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FIG. 2: (Color online) (a) DT from the numerical results compared with the experimental results in single-layer graphene on 
SiC>2 in Ref. 0- The calculation without the RI phonons (dash-dotted curve) is also plotted here. The results are normalized 
as Ref. 0- The black dotted curve represents the exponentially fitted curve of the DT in the time range of 2-4 ps. In the 
inset we plot IV calculated from Eq. (|A.4[) with (solid curve) and without (dashed curve) the second term in the equation, (b) 
Temporal evolutions of the temperatures of carriers and phonons. Here the carrier temperature is shown to be fitted well by 
To + Ta exp(— Trt) (black dotted curve) for t > 1.5 ps with Ta = 1240 K and Ft = 0.16 ps -1 . The inset zooms the time range 
4-8 ps. (c) Angular integration I\(E) from Eq. (|19p for lYo and Tto phonons as function of the carrier energy E. The two 
black dotted lines indicate E = and cur, respectively. 



be rewritten into 



Te-Tp 1 dfi+ 

r 2fc B T e 2 k B T e dt 



(17) 



by considering that T e can be fitted with T e = Tq + 
Ta exp for t > 1.5 ps, as shown in Fig. [2fb). We 
plot T r calculated from this equation with and without 
the second term in the inset of Fig. [2ja) . It is seen that 
the first term in the equation is dominant and exhibits 
a peak at T e = 2Tb = 600 K. In the time range inves- 
tigated here, T e is larger than 2Tq and thus r r shows 
a slight increase with increasing t (decreasing tempera- 
ture). Nevertheless, in the slow relaxation regime inves- 
tigated in the experiment, i.e., 2-4 ps, the exponential fit 
of DT is still acceptable because the corresponding T r 
only changes by about 8%. 

Figure [U[b) also shows that the temperatures of 
phonons in different branches are very different, originat- 
ing from the different carrier-phonon scattering strengths 
in different phonon branches. Interestingly, the temper- 
atures of the Tto and Tlo phonons differ much even 
though their carrier-phonon scattering matrices are very 
similar [sec Eq. ([5])]. This phenomenon comes from their 
different angular dependences, which is neglected in the 
simple model in the experimental works^ To show this 
more clearly, we utilize the conditions that the carrier 
distribution is isotropic and the distribution of optical 
phonons is independent of q (i.e., n q A = n\), and rewrite 
Eq. CEU) as 



dtn x \ m = — / dEI x (E)[f<(E-Lu r )f>(E)n x 
f>(E-u>r)f<{Ef)(nx + l)], (18) 



'° P iVph J-c 



where />•<(£) = /> k < 
tion I\(E) is given by 



\f l rXHH l 1 2 

l Ju k^,k'i/' 



and the angular intcgra- 



6{e Vv >-E)5(e kl/ -E+(Or). (19) 



I\{E) for Tlo and Tto phonons arc plotted as function 
of E in Fig. He). One finds that I X (E) for r LO (r TO ) 
phonons is larger than the other one in the regime E > ojy 
and E < (0 < E < Wr)- For the investigated excita- 
tion, w pu is much larger than ujy- Thus the intraband 
carrier-phonon scattering (corresponding to E > wp and 
E < 0) dominates the carrier-phonon thcrmalization. 
Consequently, the scattering strength of Tlo phonons 
is stronger and the corresponding temperature is higher. 
Another more interesting phenomenon shown in Fig.^b) 
is that the carrier temperature can be even lower than 
the hottest phonon one. This can be understood as 
follows: When the hot carriers are in equilibrium with 
the hottest phonons, the cooling of the carrier is due 
to the energy exchange with the other colder phonons; 
whereas the cooling of the hottest phonons comes from 
the anharmonic decay of hot phonons. As shown above, 
the carrier-phonon thcrmalization is faster than the hot- 
phonon decay. Thus the temperature of carrier decreases 
faster and hence becomes lower than that of the hottest 
phonons. We also discuss the contribution from the RI 
phonon, which is neglected in the literature^ by plot- 
ting the DT from the calculation without the RI phonons 
in Fig. (5Ja) . It is seen that the exclusion of this phonon 
scattering makes a marked difference, indicating that the 
RI phonons are very important to the cooling process. 

We then investigate the temporal evolutions of DT 
and temperatures of carriers and phonons in graphene 
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FIG. 3: (Color online) (a) Temporal evolution of the DT from 
the numerical results compared with the experimental data 
on SiC substrate extracted from Fig. 1 in Ref. |8j- The results 
are normalized as that paper, (b) Temporal evolutions of 
the temperatures of carriers and phonons in the two-layer 
graphene. 



on 6H-SiC. The results are compared with the experi- 
mental data reported by Wang et alM (Fig. 1 in that pa- 
per) in Fig. |3Ua)4^ As mentioned above, the RI phonons 
are important only when the number of the graphene 
layers is small. Therefore, we include the carrier-RI- 
phonon scattering for the two-layer sample (Sample C 
in Ref. d) and exclude it for the 55-layer sample (Sam- 
ple A). As presented in Ref. d, w pu = uj pl = 1.6 eV 
and the average equilibrium carrier densities for 55-layer 

and 



and two-layer samples are taken to be lxlO 11 cm 2 



6x10 cm , respectively. N e; 



8.6 x 10 1 



and 



I a = 2 mJ/m 2 are consistent with the estimated values in 
the experiment. Then the corresponding initial temper- 
atures and chemical potentials are T = 4193 K, [j$ = 
— 1061 meV and fjP_ — 1099 meV for 55-layer sample and 
T e ° = 3512 K, = -663 meV and ^° = 825 meV for 
two-layer sample. The fitting gives _E max = 0.8 eV and 
— 2.5 ps. 

Figure 3(a) shows good agreement between our calcu- 
lations and the experimental data in both samples, indi- 
cating that the contribution from the RI phonons can be 
responsible for the difference between the DTs in these 
two cases. It is noted that to fit the results, the carrier- 
optical-phonon interaction parameters (-Dp) and (D^) 
are chosen to be the same as those adopted in Ref. 8 
but twice of those obtained from the density functional 
calculations ! 15 ' 16 In fact, the values of these parameters 



are still in debate i 15 ' 16 ' 47- — Especially, the influences of 
the electron-electron correlatio n 47 ' 50 and the interlayer 
coupling 51 on the carrier-phonon interaction parameters 
in the epitaxial multilayer graphene are still unclear. 
Thus these parameters can be sample dependent. We 
also show the evolution of carrier and phonon tempera- 
tures in the two-layer graphene in Fig. G^b). From our 
results in Fig. and (b) , it can be seen that the behav- 
iors of carriers and phonons of graphene on a SiC sub- 
strate are similar to those on a SiC>2 substrate, i.e., a fast 
DT relaxation of hundreds of femtoseconds followed by a 
slower picosecond one as well as a lower carrier temper- 
ature compared with the hottest phonons. In addition, 
the DT in the slow relaxation regime decays exponen- 
tially in a large time range. This is because T e is around 
2Tq where T r varies mildly with T e as shown in the inset 
of Fig Ha). 



C. Excitation-density and probe-photon—energy 
dependences of the slow DT relaxation 

In this subsection, we first compare the calculated DT 
with the experimental data extracted from Fig. 3(b) in 
Ref.dfor different excitation densities in Fig. Efla) 4 s - Here 
the parameters are the same as those in Fig. [3] unless 
otherwise specified and the RI phonons are not included 
since the layer number of the sample is 16. The fittings 
give the excitation densities N ex = 2,5 and 8x 10 11 cm~ 2 
as well as E max = 0.6, 0.7 and 0.8 eV for the experimental 
results with pump-pulse energies being 2.5, 6 and 9.8 nJ, 
respectively. From this figure, one finds good agreement 
between our numerical results and the experimental data 
for all three excitation densities. 

By assuming that E max is a linear function of N ex 
and fitting the above values of E max and N ex , we obtain 
E max — 0.033^2, + 0.533 (E max and N ex are in units of 
eV and 10 11 cm -2 ). With this relation, one can obtain 
the evolution of DT for other excitation densities. Here 
we concentrate on the relaxation rate in the time range 
of 2-3 ps, since the evolution of DT in this regime shows 
a good exponential decay in the whole excitation den- 
sity range in this investigation. The results are plotted 
as dots in Fig. EJb). It is seen that mild valleys ap- 
pear in the N ex dependence and the excitation densities 
where valleys appear tend to be lower for smaller w pr . 
To better understand this phenomenon, we also plot the 
results from the approximate formula Eq. (IA.4I) (solid 
curves) in Fig. 0Jb). Here T e , dT e /dt and d/J,+/dt in 
Eq. (|A.4[) are chosen to be the ones at the middle of 
this time region, i.e., t = 2.5 ps for each N ex . One can 
see that the results from Eq. (|A.4[) agree very well with 
those from the kinetic equations and exhibit valleys at 
the same N ex . The scenario of these valleys is as follows. 
With the increase of excitation density, the carrier and 
phonon temperatures increase. Thus the electron-hole 
recombination and the cooling of carrier-phonon system 
both accelerate due to the enhanced carrier-phonon and 
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FIG. 4: (Color online) (a) DT from the numerical calcu- 
lations compared with the experimental data [Fig. 3(b) in 
Ref. The pump pulse energies in the experiments are 2.5, 
6 and 9.8 nj from the bottom to the top. (b) Slow relaxation 
rates of the DT from the numerical calculations (dots) and 
the approximate formula Eq. (IA.4|) (solid curves) as function 
of the excitation density N ex for different probe-photon en- 
ergies. The DTs calculated from the first term in Eq. (|A.4|) 
(dashed curves) are also plotted here. Note that the scale of 
I a corresponding to iVe^ is on top of the frame. The arrows 
show the positions of the valleys, (c) Slow relaxation rate of 
the DT as function of the probe-photon energy for different 
excitation densities. The results are fitted with linear func- 
tions (solid lines). 



phonon-phonon scatterings. As a result, —dfi + /dt and 
—dT e /dt in Eq. (|A.4[) increase with N ex . Furthermore, 
our calculations show that —d[i+/dt almost increases lin- 
early with N ex and —dT e /dt ~ N® x in the excitation 
density range investigated here. The excitation-density 
dependence of T e is more complex: for T e ~ N" x , a is 
around 0.18 for 1 x 10 11 cm~ 2 < N ex < 5 x 10 11 cm" 2 and 
then decreases slowly with increasing N ex and reaches 0.1 
for N ex = 1.3 x 10 12 cm~ 2 . Therefore, the first term in 



Eq. (|A.4p first decreases and then increases with increas- 
ing N ex , while the second term increases monotonically 
with N ex . Under the joint effects of these two terms, the 
relaxation rate shows a valley at the excitation density 
lower than that solely from the first term [dashed curves 
in Fig. HJb)]. Also by considering that the contribution 
of the first term decreases with decreasing w pr , the val- 
ley moves to lower N ex when w pr becomes smaller. We 
also present more detail about the probe-photon-energy 
dependence in Fig. |Uc). It is seen that the relaxation 
rate increases linearly with the increase of probe-photon 
energy cj pr . This can also be understood via Eq. (|A.4[) if 
one notices that /i+ and T e are independent of w pr . 



IV. CONCLUSION 

In conclusion, we have microscopically investigated 
the dynamics of nonequilibrium carriers and phonons 
in graphene by solving the kinetic equations with the 
carrier-phonon and the carrier-carrier scatterings explic- 
itly included. The hot-carrier Fermi distribution is found 
to be established within 100 fs. Furthermore, the temper- 
atures of electrons in conduction and valence bands are 
shown to be very close to each other even when the exci- 
tation density is comparable with the equilibrium carrier 
density. This is shown to be due to the strong inter- 
band Coulomb scattering. Moreover, the temporal evo- 
lutions of the DT obtained from the kinetic equations 
agree well with the experimental results-^ for different 
graphene layer numbers and excitation densities, with a 
fast relaxation about hundreds femtoseconds followed by 
a slow picosecond one presented. Based on the results 
of the evolutions of carrier and phonon temperatures, 
we find that the mechanisms leading to these two relax- 
ations are the fast carrier-phonon thermalization and the 
hot-phonon decay, respectively, which is consistent with 
the conjecture in the previous experimental works^— We 
also show that the temperatures of the hot phonons in 
various branches are very different due to their differ- 
ent carrier-phonon scattering strengths. Particularly, in 
spite of the similar carrier-phonon interaction matrices, 
the scattering strengths of the TO and LO phonons near 
the r point are very different due to their different angu- 
lar dependences. In addition, the temperature of carriers 
can be lower than that of the hottest phonons. This 
comes from the fact that the phonon temperatures are 
different for different branches and the hot-phonon de- 
cay is unimportant during the phonon thermalization. 
Our calculations also show that the contribution of the 
RI phonons is important in the relaxation process. 

Finally, we investigate the excitation-density and the 
probe-photon-energy dependences of the slow DT relax- 
ation rate. The relaxation rate is found to exhibit a mild 
valley in the excitation density dependence. This phe- 
nomenon comes from the competition among the increas- 
ing carrier temperature and the accelerating electron- 
hole recombination and carrier-phonon cooling with an 
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increase of the excitation density. We also show that Also considering that w pr is much larger than /i+ , /i_ and 

the slow relaxation rate is linear with the probe-photon ksT e , AT pr = T pr — Tp r can be expressed as 

energy. 
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Appendix: Approximate formula of slow relaxation 
rate of DT 



We present the derivation of the approximate formula 
of the slow relaxation rate of DT here. In this time re- 
gion, the Fermi distribution has been established. Thus 
the distribution functions for electrons in the conduction 
and valence bands are f c — {exp[/3(w pr /2 — + 1} _1 



and /„ = {exp[ 



r /2 - M _)] + 1}- 



l/(k B T e ). Since N lay a(uj pT ) v /// /eo/(l 



T 



Nu y e 2 V/io/ e o 
2(1 + n rcf ) 
1 



Thef 



1 



with 
) < 1. 



3 /3(w pr /2-p+) _|_ i 



(A.l) 



w e 2 iV lay VV^ e -^ pr/2(e ^ + + e _^_^ (A 2) 
2(l + ra re fj 



Thus, the relaxation rate of DT is given by 



1 



AT pr dt 



(AT pr ) 



uj^dfi _ dp+ + jg-) 

2 dt P ^ e ^++^-) + l' 



(A.3) 



In the excitation density range investigated here, 
\dfi+/dt + dfi-/dt\ <C Therefore, one has 



w pr dT e 1 dfi-\ 



2k B T} dt k n T„ dt 



(A.4) 
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